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ABSTRACT 

The motion of a charged particle in a linear electromagnetic device 
can often be analyzed by using the transport matrices. For an electron 
storage ring, this technique has been applied to yield fruitful results 
such as the trajectory of the particle distribution center and the beam 
sizes and shapes in phase space. Coupling effects among the horizontal, 
vertical and longitudinal motions are included in a straightforward manner. 
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I. INTRODUCTION 

Ignoring nonlinear perturbations, the equilibrium particle distri- 
bution, in an electron storage ring is gaussian in the phase space with 
canonical coordinates X= (x,x f ,y ,y 1 ,z , <5) , where x, y, and z are the hori- 
zontal, vertical, and longitudinal displacements relative to the center 
of particle distribution; x 1 and y' are the corresponding conjugate 
momenta; and 5 = ~ i s the relative energy deviation. ,Z,J » 4 To describe 
this particle distribution, it is necessary and sufficient to specify 
the distribution parameters which include the 6 first-moments <x > and the 21 
second-moments <x iX ^>, i,j = 1,2, ...6 with x.'s denoting the canonical 
coordinates of the state vector X. Various storage ring quantities, such 
as the closed orbit trajectory of the beam, the beam width, beam height, 
bunch length, tilting angle in the x-y plane, beam energy spread, etc., 
are directly given by those distribution parameters. Since the closed 
orbit and the beam size and shape must fulfill the requirements set by 
the design goals of storage ring, we must be able to obtain the distri- 
bution parameters for any storage ring lattice being designed. The 
matrix formalism presented in this paper will serve as a straightforward 
method of obtaining the beam distribution parameters under a wide variety 
of storage ring operating conditions. 

II . BACKGROUND 

As an electron circulates in the guiding magnetic field, it emitts 
photons due to synchrotron radiation. The energy loss of an electron 
due to emitting a photon excites the transverse and longitudinal occil- 
lations of the electron in its subsequent motion. In general, this 
effect of quantum excitations causes the beam size to diffuse in all 
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three dimensions. As the beam emittances grow due to diffusion, the 

2 3 5 

effect of radiation damping * ' becomes significant and an equilibrium 

beam distribution is reached as the quantum diffusion and radiation 

damping balance each other. 

Assuming an . ideal storage ring and ignoring any coupling between 

x, y and z dimensions, most beam distribution parameters simply vanish 

and the remaining parameters can be obtained by relatively elementary 

2 3 

considerations ' . However, in practice, some coupling effects may not 
be negligible. For example, x- and z-motions are coupled if the rf cavi- 
ties are located at positions with finite energy dispersion; x- and y- 
motions are coupled if there are skew quadrupole fields in the storage 
ring. One possible way to include coupling in the evaluation of distri- 
bution parameters is to assume that the coupling effects can be approxi- 
mately described by a set of coupling coefficients which specify the 

4 

coupling strength averaged over one revolution of the storage ring. Distri- 
bution parameters are then obtained by solving the corresponding Fokker- 

1 2 

Planck diffusion equation. * Another method ignores couplings to the 
longitudinal z dimension and the relevant distribution parameters are ob- 
tained by a careful analysis of the balance between quantum diffusion and 

7 8 

radiation damping. ' In the above theories, it has been necessary to 

9 

introduce the horizontal and vertical energy dispersion functions in 
order to carry out the analysis . 

In the following we describe a general framework which determines 
the distribution parameters in a coherent manner and with the above- 
mentioned coupling effects fully taken into account. In this method, each 
linear element in the storage ring lattice is represented by a 6x6 
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TRANSPORT matrix 1 ^ which transforms the state vector X as an electron 
passes through the element. Nonlinear elements are approximated by 
linearizing their electromagnetic fields around the center of particle 
distribution. Knowing the TRANSPORT matrix transformations around the 
storage ring, the distribution parameters can be obtained from the eigen- 
values and eigenvectors of some proper matrices yet to be described. In 
particular, no knowledge of the energy dispersion functions is needed. 

III. TRAJECTORY OF THE BEAM CENTER 

Relative to the design orbit, the motion of an electron can be de- 
composed into two components, and X, with Xq describing the trajec- 
tory of the beam distribution center and X being the oscillatory devia- 
tion from Xq. The fact that the trajectory of the beam distribution 
center does not coincide with the design orbit (i.e. (^) could be due 

to several causes: (i) A quadrupole magnet misaligned in its transverse 
position is the same as superimposing a dipole magnetic field to the 
magnet. This accidental dipole field kicks the beam orbit transversely; 
(ii) Additional dipole magnets may exist in the lattice for orbit correc- 
tion or beam injection purposes; (iii) The energy loss of an electron 
due to synchrotron radiation is compensated by the acceleration provided 
at the rf cavities. Since the rf cavities are usually not located at 

places where synchrotron radiation occurs, the beam distribution center 

1112 

can deviate from the, design trajectory as a result. » 

Existing methods of obtaining Xq treat each orbit distortion 
mechanism individually and ignore the possibility of couplings among x, 
y, and z motions. In the following of this section, we present a method 
which includes various orbit distortion mechanisms simultaneously and 
taking into account the coupling effects. 
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We first form a 7 -dimensional vector V by adding to a seventh 

13 

component which is always given by unity. The vector at the exit of 
a given lattice element is linearly related to the vector at the 
entrance by V f = MV^ The 7x7 transformation matrix M can be obtained 
from the electromagnetic fields in the lattice element. For simplicity, 
we assume that all lattice elements aside from the drift spaces are short 
compared with the focal lengths so that a thin-lens approximation applies. 
Generalization to the thick-lens case is straightforward. Under this condi- 
tion, the diagonal elements of M are always equal to unity and off -diagonal 
elements except those listed in Table 1 vanish. In Table 1, we have de- 
fined £=length of the lattice element, Bp = particle rigidity, 
2irR = circumference of the storage ring, = design particle energy, 



C = 
Y 



4tt 



T~ r 0 / ( m o c ) with r the classical electron radius and m the elec- 
o e e e e 

3 " 

tron mass; and for the rf cavities, V = peak applied voltage, h = har- 
monic number, <j> = synchronous phase. The upper left 6*6 corners of 
these matrices are the thin-lens version of the usual TRANSPORT matrices. 10 

The trajectory of the particle distribution center, described by 
V(s) = (X Q (s), 1) = (x Q , x^, y Q , y^, z Q , 6 Q , 1), with s indicating the 
position around the storage ring, must satisfy 

X 0' X 0' Y 0' y 0 ( s+2ttR ) = x o» X 0» y 0' y 0 (s) ' 

net energy change per turn =0 (1) 

and 

net path length travelled per turn = exactly 2irR, 
where the last condition is imposed by the fact that the rf cavity frequency 
is accurately fixed. To find V(s), we first look for the transformation 
matrix for one revolution, W(s) , by multiplying the M-matrices of all 
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lattice elements successively from s to s + 2?rR. V(s) is then simply 
given by the eigenvector of W(s) with eigenvalue 1, i.e., 

W(s) V(s) = V(s) . (2) 
That one of the eigenvalues of W(s) is equal to 1 can be easily proved. 
The trajectory of the particle distribution center obtained here, X^s) , 
will also be referred to as the closed orbit of the beam. 

Effects on closed orbit due to nonlinear elements . can only be treated 
by iteration. A sextupole magnet is thus represented by a horizontal and 
a vertical kicker with the nontrivial transformation matrix elements 



SEX: m 26 = - m 51 = -m^ = y <*g- y g) , (3) 

~ m 46 = m 53 = m 47 = A x 0 y 0 ' 

where A = ~~ 3 2 B /3x 2 is the strength of the sextupole field. Nonlinear 
Bp yj 

effects due to radiation energy losses in quadrupole, skew quadrupole and 
sextupole magnets can also be included by inserting the following addi- 
tional matrix elements: 



SEX: m 6? = -C E^A 2 (xg + yg) 2 /8-rr^ 



y 

SKQ: ^ - -C t E3^ ^) (xg + y§)/2«£ (4) 

QUA: m 67 = -c^A^j (jt 2 + y2)/2^ . 

The iteration steps are: (i) calculate V without nonlinearities , 
(ii) evaluate transformation matrices for nonlinear elements according 
to Eqs. (3) and (4), (iii) Recalculate V with the new matrices, (iv) iter- 
ate. A few iterations should be enough for convergence. The total energy 
loss per turn due to synchrotron radiation can be obtained from 
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U_ = 52 eV (sin<j> + 77 z n cos<}> 1 
0 \ s R 0 s/ 

CAV 



IV. RADIATION DAMPING CONSTANTS 



2,3,5,7,8 



(5) 



Once the closed orbit is obtained, a particle's motion can be de- 
scribed by the state vector X=(x, x', y, y T , z, S), which is the deviation 
from closed orbit. The transformation matrices for X are given by the 
upper left 6x6 corner for lattice elements listed in Table 1. The trans- 
formation for sextupoles, linearized around the closed orbit, is a combi- 
nation of HK, VK, SKQ and QUA; 
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(6) 



From here on, we will concentrate on these 6x6 TRANSPORT transformations. 
By successively multiplying all M-matrices from s to s+2-rrR, one gets the 
transformation for one revolution, T(s), which will be extensively used 
in the next section. To obtain the radiation damping constants, however, 



further modifications on the matrix elements are necessary. 
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As an electron passes through an rf cavity, it receives a boost in its 
longitudinal momentum, while keeping the transverse momentum unchanged. 



Consequently, the slopes x' and y f of the electron T s trajectory is reduced 

u 
E, 



in magnitude by a factor of 1 - , where u = eV ^sin<j) g + ^ z^ cos^jis 



the energy gain of the electron at the rf cavity. This effect provides 
the transverse radiation damping to the motion of an electron. In the 
matrix formalism, it is taken into account by replacing the following 
diagonal element in the CAV matrices : 

m 22 = m 4A = 1 ~ ( Sin *s + I z 0 COS *s) ■ (7) 

The longitudinal radiation damping comes from the fact that an 
electron with higher energy tends to lose energy by radiating more in a 

magnetic field. This results in the following changes: 

2 

HB, VB, HK, VK: ^ = 1 - C^^J hi - 
QUA, SKQ: = 1 - A q (xg + y§) ■ (8) 

SEX: m 66 = 1 - \ A g (x§ + y§) , 
/ 3B \ 2 

where A = C E^ ~- hi and A = C E^A 2 ^ + y 2 )/2iT£. In addition, 

q y 0\Bp 9x,y/ s y u U J U 

energy radiated in quadrupole, skew quadrupole and sextupole magnets also 

depends on the transverse displacements of an electron. This contributes 

to the mixing among the transverse and longitudinal damping effects. 

Upon linearization, it yields 



QUA, SKQ: m^ = -A q x Q , ^ « -A q y Q 
SEX: m 61 = -Ax Q , m 63 = -A^ # 



(9) 
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After making the modifications (7), (8) and (9), one calculates once 
again the transformation matrix for one revolution, D(s), by multiplying 
M-matrTces from s to s+2ttR. The eigenvalues of D(s) are expressed as 
exp (-a^±i2Trv^) , k=I, II, III, where k is the mode index, is the tune 
of mode k and o^'s are the radiation damping constants. In particular, 
if any one of the ct^ T s is negative, the motion will be unstable. By eval- 
uating the eigenvalues of D(s), we have thus derived a method of calcu- 
lating the radiation damping constants in the presence of various cou- 
pling effects. Furthermore, it follows from the property of eigenvalues 
that 

-2(a I + a 1T + a TII ) = det ( D ( s )) . (i 0 ) 

e 

ii 5 
For weak damping with |ct, | << 1, it reduces to the usual sum rule 

a I + a II + a III = ^0^0 * " (ii) 

In the special case of no closed orbit distortion and no coupling, the 
damping constants are given by 

a x =a y =U 0 /2E 0 , ^ = U Q /E 0 . (12) 

V. BEAM SIZES AND SHAPES 
Our method can be used to determine quantities like the horizontal 
and vertical beam sizes, the tilt angle in the x-y plane, the natural 
bunch length and energy spread, etc. To do this, let us consider a pho- 
ton being emitted at position s^ with energy deviating from the mean value 
by a random amount SE. Let T(sq) be the coupled 6*6 transformation matrix 
for one revolution obtained in the previous section without radiation 
damping. The eigenvalues, A^, and eigenvectors, E k (s Q ), of T(s Q ) with 
k = ±1, ±11, ±111 are defined by 



T(s o ) V s o> - \W 



E k< 8 0> = E -k (s o ) 



(13) 



±i2TTV, 



A = e 

±k - 
Following the photon emission, the subsequent motion of the electron is 



described by 



X(s) =2 A^Cs) 
k 



> s - s 0 



(14) 



where E^_(s) is the eigenvector of T(s) obtained from E^s^) by the matrix 
transformation from Sq to s. Equation (14) satisfies the initial condi- 
t ion 



0 
0 
0 
0 
0 



L" 6E/E 0- 



= I Vk <8 0> 

k 



(15) 



where the left-hand side is the impulse perturbation to electron state 
due to the photon emission event. It should be pointed out that the 
transverse dimensions of the electron are not excited at the instance 
of photon emission and consequently no knowledge of the energy dispersion 
functions is needed. 

9 

From the symplecticity condition of T(s), i.e. , 



TST = S 



(16) 



where a tilde means taking the transpose of a matrix and 
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one can prove that 



E.SE. = 0 unless j = -i 



We will normalize the eigenvectors so that 



(17) 



E k SE k = i , k = I, II, III . (18) 

This normalization condition is preserved as a function of s due to the 
symplecticity of T(s). Using Eqs. (17) and (18), Eq. (15) yields 



. . 6E , . 

\ = - 1 1^ E k5 (S 0 } 



(19) 



th 



where E^ means the i component of the vector E^. Assuming all photon 
emission events are uncorrelated , one obtains the quantum diffusion rate 
of JA^I 2 by averaging Eq. (19) around the storage ring: 



dt 



|A k |2>=^/ ds <N|^>| Ek5 (s)|2 . (20) 

2 3 



In Eq. (20), N is the number of photons emitted per unit time and 

<N |V = 2C tY 5/ |p(a)|3 , 



(21) 



J 0 



where C = 

Li 



55 



r -fT/m with 1i the reduced Planck's constant, y is the 
48/3 6 e 



relativistic factor and p(s) is the bending radius. 
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So far we have ignored the radiation damping which, when taken into 
account, gives an additional contribution 

iHV 2 * ■ -^f <|A K |2> • (22) 

with Tq the revolution .time and a^. the radiation damping constants found 
in the previous section!.. The equilibrium values of <)a^| 2 > is given by 
a balance between quantum diffusion and radiation damping, which gives 

<l\| 2 > = <|A_ k l 2 > = C L ^/ds-^^ 2 • (23) 

k p(s) P 



It follows from Eq. (14) that the particle distribution parameters at 
position s are given by 



< x iV (s) = 2 2 ^ 2> Re 

k=l,ll,Hl 



E ki (s) E£ (s) 



(24) 



Equations (23) and (24) are our final expressions. . The tilt angle 6 of 
the x-y beam profile relative to the horizontal axis can be found from 

tan26 = Jt% 2> (25) 
and the transverse beam area (for luminosity calculation) is given by 

A = -rr V<x 2 xy 2 >-<xy> 2 . (26) 
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